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NUMERICAL SIMULATION OF NOISE PROPAGATION IN JET ENGINE DUCTS 


by Kenneth J. Baumeister and Edward C. Bittner 
Lewis Research Center 

SUMMARY 

A finite difference formulation is presented which should be useful in the study of 
acoustically treated inlet and exhaust ducts used in turbofan engines. The difference 
formulation can readily handle acoustic flow field complications, such as axial varia- 
tions in wall impedance and cross-sectional area, that would occur in a sonic inlet. 

The formulation is both mathematically rigorous and convenient to use . 

In formulating the difference solutions, the continuous acoustic field is lumped into 
a series of grid points spread uniformly throughout the acoustic field. At each grid 
point in the acoustic field, the pressure is separated into its real and imaginary terms. 
Then, the two-dimensional Helmholtz equations for both the real and imaginary pres- 
sures at each grid point are expressed in difference form. The solution is obtained by 
coupling the difference equations for the real and imaginary pressures with the impe- 
dance conditions at the boundaries. The solution yields the two-dimensional distribution 
of pressure, velocity, and intensity level throughout the duct. 

The coefficient matrix of the acoustic difference equations is presented and dis- 
cussed. Then, example solutions are presented for the propagation of sound in a one- 
dimensional straight hard -wall duct and in a two-dimensional straight soft-wall duct 
without steady flow. These numerical solutions showed favorable comparison with the 
exact analytical results. The examples are used to establish rules of thumb for choosing 
grid sizes for various frequencies and duct lengths. 

Because the solution matrix for the acoustic flow field is not positive definite, con- 
ventional iteration techniques cannot be used to solve the difference equations. Before 
the finite difference formulation can be applied to more complicated problems which re- 
quire a large number of grid points, such as the sonic inlet, iteration techniques will 
have to be developed to overcome present grid -size limitations. 



INTRODUCTION 


The level of noise currently being proposed for new aircraft requires the use of 
acoustic suppression in the inlet and outlet ducts of their engines. At the present time, 
there is a need for more flexible suppressor design techniques which can handle acous- 
tic flow field complications in engine ducts, such as axial variations in wall impedance 
and cross-sectional area, that would occur in a sonic inlet. To meet this need, the 
present report develops a numerical finite difference technique which can be used in the 
prediction of sound attenuation in turbojet engine ducts. In addition, acoustic wave 
propagation is important in ultrasonics and underwater technology, as well as in con- 
ventional noise abatement problems in architectural engineering. The techniques pre- 
sented herein might also be applied to these other areas of acoustics. 

Morse (ref. 1), Cremer (ref. 2), and Rice (ref. 3) have presented analytical solu- 
tions for the transmission of sound in ducts of infinite length, uniform cross section, 
constant wall impedance, and no flow. The theoretical solutions were later extended by 
Rice (ref. 4), Eversman (ref. 5), and Lambert and Tack (ref. 6) to acoustic propagation 
with uniform subsonic flows. Pridmore -Brown (ref. 7), as well as other authors (refs. 

8 to 17), has shown that shear flows (boundary -layer effects) can have a significant in- 
fluence on the attenuation of sound in ducts. All these works use the technique of separa- 
tion of variables to obtain a complex eigenvalue formulation in which the eigenvalues are 
determined by various numerical techniques. 

The present report formulates a solution for the propagation of sound in ducts with 
arbitrary wall impedance by a finite difference technique. A search of the literature 
failed to uncover any previous work in this area. The difference solutions bypass the 
conventional eigenvalue problem with its associated modes. Asa result of the difference 
formulation, the propagation of noise is treated as a diffusion process analogous to prob- 
lems in thermodynamics involving heat flow or those in fluid dynamics involving the 
transport of vorticity. 

Immediately following the mathematical development of the difference technique, 
two examples are presented to illustrate the numerical technique. These examples have 
analytical solutions; therefore, a direct comparison between the numerical and analyti- 
cal results can be made. These examples are also used to establish rules of thumb for 
choosing grid sizes for various frequencies and duct lengths. By choosing the maximum 
acceptable grid size, the amount of computer time necessary to obtain a solution to the 
finite difference equations can be minimized. 

The first example presents a solution for the sound propagation in a one -dimensional 
hard -wall duct with the impedance specified at the duct exit. The second example treats 
the propagation in a two-dimensional straight soft-wall duct without steady flow. Also, 
the limitations of the present theory are discussed and suggestions for future work are 
presented . 
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GOVERNING EQUATION AND BOUNDARY CONDITIONS 


The governing differential equation and boundary conditions are introduced in this 
section of the report. First, the equations are presented in dimensional form. Next, 
the equations are made nondimensional, and the significance of the dimensionless pa- 
rameters is briefly discussed. 

The basic equations for the scalar pressure field and the velocity field are pre- 
sented, as well as the appropriate expressions for the wall impedance condition and 
acoustical intensity . 


Scalar Pressure Field 


The continuity, momentum, energy, and state equations of a perfect gas in a duct 
reduce, in the case of no steady flow, to the linear two-dimensional wave equation 


V' 2 P' - 1 a 2 p» 


2 2 i.2 

c 3t 


(1) 


In equation (1), c is the speed of sound, t is time, and P’ is the pressure fluctua- 
tion which is assumed to be small compared to the ambient pressure. The prime (’) is 
used to denote a dimensional quantity. (These and all other symbols used in this report 
are defined in appendix A . ) 

The assumptions involved in the derivation of equation (1) are given in most acous- 
tic texts and are not discussed herein. For ease in illustrating the numerical tech- 
niques, we have chosen this simple form of the wave equation rather than the more com- 
plete form, which contains Mach number and shear flow effects such as given by Mungur 
and Gladwell (ref. 8). 

We now make a further simplification by working in rectangular coordinates. In 
this case, equation (1) can be expressed as 


g> 2 P’ d 2 P’ _ J_ 3 2 P’ 
Sx’ 2 3y’ 2 c 2 9t 2 


( 2 ) 


where 


P’ - P’(x’,y\t) 


(3) 
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This coordinate system is illustrated in figure 1 . 

We are concerned primarily with the steady -state solutions of the wave equation, 
equation (2); consequently, we neglect the transient startup effects in the system. Thus, 
for steady state we assume a solution of the form 


P*(x' , y' > t) = p'(x',y')e +ia,t 


(4) 


where in general 


p’(x',y') = p'^(x',y') + ip’^(x',y’) (5) 

The (1) and (2) superscripts designate the real and imaginary pressure components. 

Both p'^ and p'^ are real -valued functions by definition. Complex quantities are 
used because of the ease of manipulation; however, only the real part of the solution will 
have physical significance. Substituting equation (4) into the scalar wave equation 
(eq. (2)) yields 

si: + 2V +k V = o (« 

ax * 2 ay ’ 2 

where the wave number k is given by 

k = — (7) 

c 

Equation (6) is commonly called the Helmholtz equation. This equation and its appro- 
priate boundary conditions will be solved numerically by the finite difference technique . 

In nondimensionalizing equation (6) we choose to write equation (6) in such a form 
that the numerical results will be compatible with present NASA Lewis acoustical design 
procedures. To do this, we define the following set of dimensionless variables as sug- 
gested by the results of Rice (ref. 4): 


y 



( 8 ) 


x 



(9) 
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p = SL 

pa 


(10) 


where H equals the height of the rectangular duct and will usually be chosen as the 
amplitude of the imposed pressure wave at x = 0. 

Substituting equations (8) to (10) into equation (6) and rearranging terms gives 


i?P + ^E+ (27TTj) 2 p = 0 (11) 

ax 2 0y 2 


where 


JH w = Hf = H 

277 C C X 


( 12 ) 


The frequency parameter rj represents the ratio of duct height to acoustic wavelength. 
The dimensionless height y ranges from zero to 1, while the dimensionless length x 
ranges between 

0< x < — (13) 

H 

where L is the length of the duct. Thus, the important dimenionless parameters de- 
scribing the acoustic attenuation are r\ and L/H. 

Equation (11) now can be broken into its real and imaginary parts. Substituting 
equation (5) into equation (11) yields 




(14) 


3x 


ay 


2 ( 2 ) 2 ( 2 ) 2 
2-E + i-2 + ( 2777 ?) V 2) - 0 


(15) 


ax' 


ay 


Acoustic Particle Velocity 

The boundary conditions in acoustics are generally given in terms of impedances 
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which relate the pressure and velocity fields at the boundary. The velocity field can be 
expressed in terms of pressure from the equation 



which for a harmonic solution of the form given by equation (4) reduces to 


u' = — V’p’ 

(17) 

pw 


In terms of the scalar velocities, equation (17) becomes 


u’ = J_iE! 

(18) 

x po> 3x' 


u* =-L 

(19) 

y ay’ 



In terms of the dimensionless parameters, equations (18) and (19) become 

°x 

where 


ap 

(20) 

ax 


ap 

(21) 

3y 


H 

- u' 

(22) 


pa 


Impedance Boundary Condition 


In the transverse direction, the acoustic impedance at the walls is defined as 

Z, = -SL at y = 0 and 1 

ii » 


( 23 ) 


6 



The impedance Z t is now replaced by the impedance ratio which is defined as 

Z. 

= — (24) 

pc 

Substituting equation (19) into (23) and introducing the dimensionless parameters yield 
the following expression for the impedance at the walls: 


= -i27T?7 - 2 - at y = 0 and 1 (25) 

3p 

3y 

or, the pressure gradients at the walls become 

3p = p a t y = o and 1 (26) 

3y ? t 


Because of convenience, we shall use equation (26) to apply the boundary constraint 
rather than equation (25) . 

It is convenient to express the reciprocal of the impedance ratio in terms of the 
acoustic conductance ratio k t and the acoustic susceptance ratio o^., that is, 



i<7 + 


(27) 


where 1/^ is called the acoustic admittance ratio. Therefore, equation (26) becomes 


Q = -2ttt7(ct + ifc J P (28) 

ay 1 


Finally, expressing p in terms of its real and imaginary terms, and equating these 
terms gives 


3p^ 

9y 


-27777 


p (1 V 




at y = 0 and 1 


ap^ 

9y 


— 2 7777 


p^ 2 ^a, + p 


(D. 


at y = 0 and 1 


(29) 


(30) 
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The impedance at the exit plane, x = L/H, leads to a similar expression 


3P 


(1) 


Sx 


= -2n7] 


x=L/H 


P^°e " 


( 31 ) 


ap 


( 2 ) 


dx 


= -2vt) 


-ix-L/H 


p^Og + p^^/f e 


(32) 


We can see that the pressures p^ and p^ are coupled by the impedance condi- 
tions represented by equations (29) to (32) . 


Entrance Conditions 


At the entrance, the dimensionless pressure p(0,y) is assumed to be of a general 


form 


P(0,y) = p (1) (0,y) + ip( 2) (0,y) 


(33) 


Generally, some real pressure profile f(y) is assumed at the entrance, that is 


p(0,y) = f(y) 


(34) 


where f(y) is normalized such that 



(35) 


The function f(y) can be chosen to map any assumed pressure profile. For the prob- 
lems treated in this report, the real and imaginary pressures at x = 0 are 


P^(0,y) = f(y) (36) 

P (2) (0,y)=0 (37) 
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For a uniform profile, as is used in the examples considered later, equation (36) be- 
comes 


P (1) (0,y) = 1 


( 38 ) 


Axial Acoustic Power 

The sound power which leaves a duct and reaches the far field is related to the axial 
intensity. The axial intensity can be expressed as (ref. 4) 


r = ^Re {p*’^ 


(39) 


In terms of the dimensionless parameters, equation (39) can be written in dimensionless 
form as 


I = — Re fp*uJ 
1 XJ 


(40) 


where 


I = p 
2 

pa 


(41) 


For a hard-wall duct with a pc exit impedance, I is identical to 1 for all frequen- 
cies and duct lengths. This result is shown in appendix B. 

In terms of the complex representation, the expression for I becomes 


I = 


,(2) 0p 


( 1 ) 


27TTJ 


0X 


" P 


(1) 0 P 


( 2 ) 


0X 


(42) 


The total dimensionless acoustic power is the integral of the intensity across the test 

i 

section 



I 

(43) 
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By definition the decrease in decibels of the acoustic power from x = 0 to x can 
be written as 


AdB = 10 


E. 


log 


10 


(44) 


FINITE DIFFERENCE FORMULATION 

The continuous system will now be reduced to an equivalent lumped -parameter sys- 
tem by means of the finite difference approximations. Instead of a continuous solution 
for the pressure, we shall find the pressure at isolated grid points, as shown in figure 2. 

The previously discussed governing differential equations can be approximated in 
difference form (ref. 18) by using one of the following three methods: (1) a Taylor ser- 
ies expansion, (2) a variational formulation, or (3) an integral formulation. For prob- 
lems where the gradient is specified along a boundary, such as in our acoustic problem, 
the integration method for generating the finite difference approximations is most con- 
venient and is used in this report. 

The difference equations for the acoustic field are presented without derivation. 

The detailed derivations of these equations are given in appendixes C and D. 


Helmholtz Equation 


The Helmholtz equation, (14), in finite difference form becomes 


-pW . - pW . 

*+-1,3 *1+1 , 3 


Ax 

Ay 


p U) - i + p^- 1 

*1,3-1 *i,J+l 


2 + 2 


(if- 


(27777 Ax)' 


P^ 1 } = 0 

*+,] 


(45) 


-p( 2 V • - p® - 

*1-1, J 1+1,3 


Axr 

Ay/ 


IP^ 2) . 1 + P< 2) - , 

ri, 3-1 *1,3+1 


2 + 2 


Ax\ 2 _ 

Ay) 


(27777 Ax)' 


= 0 
1,3 


(46) 


where 


p i,j - p i!i + ip ‘!i 


( 47 ) 


and where i corresponds to x and j corresponds to y. Note, the frequency term 
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o 

(27ttj Ax) subtracts from the term which will occupy the main diagonal element of the 
coefficient matrix; consequently, for sufficiently high frequencies or spacing param- 
eters, the matrix will no longer be positive definite (ref. 18, p. 68). Asa result, con- 
ventional iteration techniques cannot be used . 


Impedance Boundary Conditions 


The impedance boundary conditions of the upper wall with j = m (eqs. (29) and (30)) 
become 


i, m 


/ \2 

1 + \Ax) + ^ ^ 


i 

i (27TT7 Ay)' 
m 2 


_ p(D _ 1 /Vf D U) 
i » m " 1 2\Ax/ i_1,m 




>< 2 > 

i,m 


1 + 


(£& + 2mj Ay (s.) JJmMl 
\Ax/ ' x 'i,m 2 


(2) _ 1/Ay 

2\Ax 


■y«. 


m 




Similar equations exist at the lower boundary (see appendix D) . The subscripts on a 
and k indicate that these quantities can vary along the length of the duct. 

At the duct exit, the impedance condition 


k = 1 
e 


(50) 


°e = 0 


(51) 


is used so that a finite duct can approximate a duct of infinite length (see appendixes B 
and E for detailed discussion). For this case, the finite difference approximations to 
equations (31) and (32) become at the exit (i = n) 
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+ 2nr\ Axp 



(53) 


Finally, special consideration must be given to the corners at the exit of the duct, 
i = n and j = m, since both the exit and side wall impedance values influence this point. 
Asa result, the impedance condition at the corner is 



3 ( 2 ) 

n, m 


1 + 



+ 2ur} Ay 



(27777 Ay) 2 
2 



( 2 ) 

n-1 , m 


- P 


( 2 ) 

n, m- 


1 



(55) 


Axial Intensity 

In terms of the difference notation, the axial intensity as given by equation (42) can 
be expressed as 
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1 


(56) 


I. . 


i] 


p.® 

2 7T7J AX [_ 1? J 





The total intensity across the test section, as given by equation (43) is written in differ- 
ence notation as 



(57) 


We shall now apply the difference equations to the problem of noise attenuation in a 
duct. 


MATRIX SOLUTION 

The collection of the various difference equations at each grid point forms a set of 
simultaneous equations which can be expressed in matrix notation as 

A • P = F (58) 

where A is the coefficient matrix, P is the pressure matrix containing the unknown 
pressures, and F is the known matrix containing the various initial conditions, where 
A, P, and F are complex in general. 

In considering solutions to equation (58), it is convenient to express this equation in 
terms of all real quantities. To accomplish this, the matrix P is written as a column 
vector in terms of the p^ and p® pressures and subdivided as follows: 


1 

0 

1 

i-H 
< 
l , 


1 

UL 

H-* 


p(l) 

c ; a, 

i 


p(2) 

ii 

F<2) 


(59) 


The upper left A 1 matrix is associated with the p^ ^ terms and the lower right A 1 

1 ( 2 ) 1 
matrix is associated with the p v ’ terms. Both matrices are identical because of sym- 
metry which exists in the governing equations. The A^ matrix has a form typical of 
those matrices found in two-dimensional heat conduction problems or neutron diffusion 
problems. The A 1 matrix has one main diagonal element and four other diagonal ele- 
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merits which represent the coupling of the central point to the four adjacent points. 

The C matrix represents the coupling that occurs between the p^ and p^ 
pressures through the impedance boundary conditions. The lower left matrix differs 
from the upper right matrix only by a change of sign. The C matrix is a sparse matrix 
with only one main diagonal term. 


To better illustrate the detailed structure of the matrix given by equation (59), con- 
sider the simple example shown in figure 3. For this simple case the detailed matrix 
structure becomes 



where the values of the coefficients a and c in this matrix are given in appendix D, 
along with the derivation of the other elements in this matrix equation. 

As shown in appendix D, the constants a and c depend on the impedance and the 
grid size. The double -subscript notation, illustrated in figure 2, has now been replaced 
by a single -subscript notation. However, the double -subscript notation will still be used 
to set up the basic form of the matrix coefficients. 
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Matrices of the form of equation (59) can be solved by elimination techniques. In 
particular, the Gauss elimination technique was used to find a solution in the example 
problems which follow. 


EXAMPLE PROBLEMS 

Two examples are presented to illustrate the numerical technique. Both examples 
have analytical solutions; therefore, a direct comparison between the numerical and ana- 
lytical techniques can be made. As mentioned in the INTRODUCTION, these examples 
are used to establish rules of thumb for choosing grid sizes for various frequencies and 
duct lengths. The first example presents a solution for a one -dimensional hard-wall 
duct with the impedance specified at the duct exit. The second example presents a solu- 
tion for the propagation of sound in a two-dimensional straight soft-wall duct. 


Example 1 - One-Dimensional Hard-Wall Duct of Infinite Length 


As our first example of the difference technique, we shall consider the problem of 
an infinite hard -wall duct with a plane pressure wave at the entrance of the duct. 

In this case, "hard wall" implies that 


V" 

O 

II 

and 

y = 1 

(61) 

ri- 

ll 

O 

y =0 

and 

y = 1 

(62) 

II 

O 

y = 0 

and 

y = 1 

(63) 


For this special case of a plane wave entrance condition, no pressure gradients will exit 
in the transverse direction in the duct; consequently, the Helmholtz representation of the 
wave equation, equation (11), reduces to the one -dimensional form 


3x 


£ + (27777) ^p = 0 
2 


(64) 


The difference form of the general equation presented earlier can be simplified by 
noting that 
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In this case, the difference form of the wave equation, equation (45), reduces to 

-pRi , ] - p,>i , j + [ 2 - < 2 ”’' pJ! ] - ° <66) 

( 2 ) 

with an identical expression for p v ' component of pressure. 

The word ’’infinite" used in the title of this example problem implies that the duct 
is sufficiently long so that no reflective waves occur in the duct. In the numerical 
formulation of the infinite duct, we can simulate the infinite duct by assuming a finite 
length of duct and applying an impedance of value pc at the exit plane of the duct (see 
appendix B for proof). As a result, k = 1 (eq. (B12)) and = 0. Thus, the exit im- 
pedance boundary conditions corresponding to equations (52) and (53) become 


Pn.’j [ l - - 2irfj Axp®. = o 

(67) 

p (2) L _ (27rr L Ax) 2 l _ J2) + 2jnj Ax p ( 2 ) - o 

ft > J 2 n-i , j ft > J 

(68) 

The entrance conditions are simply 


p (1) (0) = 1 

(69) 

o 

n 

CD 

04 

& 

(70) 


In solving the resulting set of matrix equations we will simplify the structure of the 
solution matrix. This simplification should be useful in future work if block iteration 
techniques (ref. 18) are applied to the solution of the two-dimensional matrix equations. 

Consider the simple case shown in figure 4, where four points are used to describe 
the pressure distribution in a finite hard -wall duct with pc exit impedance. In this 
case, the general matrix equation, equation (59), reduces to the following form: 




where 

a 1 = 2 - ( 27 ttj Ax) 2 

a = 1 _ (2ih7 Ax) 2 
2 2 

= 2nri Ax 

As seen in equation (71), the basic structure of the solution matrix is similar to the 
structure shown in equation (60). The A^ matrices become tridiagonal and the coupling 
C matrix has only one element. 

The detailed structure of equations (60) and (71) are arbitrary, since we are free to 
construct the column vector P in a variety of ways. Therefore, for one -dimensional 
problems, we shall now redefine the column vector P to be of the form 
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( 2 ) 

where the p v ' components have been inverted. For the special case of n = 4, as 
shown in figure 4, , 
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Expressing the general matrix equation (58) in terms of this new column vector gives 
for our simple example the following matrix equation: 


a i ■’ ! 


■pi 1T 


1 

- 1 - 1 ! 




0 

, 0 





i 

►— » 

A> 

1 

h-* 

/ 


p) 1 ' 


0 

1 

-1 a 2 , -c 




0 


X 


_ 


c | a 2 -1 

1 


(2) 

?4 


0 

1 

1 - 1 -1 
\ 1 


P (2) 

P 3 


0 

N 1 -1 a x -1 

1 


p?> 


0 

1 -1 a 

L i 1 a jJ 




0 


where aj, a 2 , and c are the same as defined in equation (71). As seen in equation (74), 
the resulting coefficient matrix is tridiagonal, which is extremely easy to solve. 

Equation (74) was applied to a duct with a height-length ratio L/H of 1 and a dimen- 
sionless frequency rj of 1. For a value of n = 12, the numerical as well as the ana- 
lytical results are shown in figure 5. As seen in figure 5, the numerical and analytical 
results are in good agreement. 

In determining the attenuation of a particular liner, we are interested in how well 
the difference formulation predicts the intensity at the duct exit. Figure 6 displays the 
numerical results for axial intensity at the duct exit for a variety of dimensionless fre- 
quencies rj. Recall that the intensity was nondimensionalized such that it is identical to 
1 for a hard -wall duct. For a sparse number of grid points n, the numerical results 
fall below the analytical value of 1. However, as the number of grid points is increased, 
the numerical value of the intensity asymptotically approaches the theoretical value of 1. 

From the numerical results shown in figure 6, we can establish the following rule of 
thumb for one-dimensional problems: 

n = 1277^=12- for Z. = «>,Z = p c (75) 

H X 1 e 
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or 


Ax = — for Z, = 00 , Z = pc (76) 

12 1 

That is, for a given L/H or 77, if we choose n greater than or equal to the value sug- 
gested by equation (75), the calculated value of I will be within 4 percent of the true 
value . 

For the two-dimensional problem with finite Z^, equation (75) will not hold. Be- 
cause the intensity can vary inside the duct in the transverse direction in the two- 
dimensional problems, more grid points will be necessary to assess the intensity in 
two-dimensional problems. Consequently, for the general two-dimensional case 

n>12 nk for Z. < °°, Z = pc (77) 

H 1 e 

Example 2 - Propagation of Sound in a Two-Dimensional Soft -Wall Duct 

In our second example of the difference technique, we shall consider the problem of 
calculating the maximum attenuation possible in a two-dimensional duct with L/H of 
0.5 and dimensionless frequencies r\ of 0. 6 and 1. We can compare the numerical re- 
sults to the analytical results found by using the techniques presented by Rice (ref. 4). 

The first step in determining the maximum attenuation in the duct is to calculate the 
attenuation at various Z -values in the impedance plane. This is shown in figure 7, where 

t=— =6 + ix (78) 

pc 

defines the impedance plane. As seen in figure 7, the peak occurs at negative values of 
wall reactance. The sound attenuation values shown in figure 7 were obtained by choos- 
ing discrete values of 9 and x throughout the 0-x plane, calculating the sound atten- 
uation at each point, and interpolating between the points to obtain smooth contours. 

The impedance associated with the peak attenuation is called the optimum impedance. 

As was previously mentioned in this report, the solution matrix for the acoustic 
flow field is not positive definite; consequently, conventional iteration techniques cannot 
be used to solve the difference equations. This is a serious limitation, since iteration 
techniques generally are quicker and require a much smaller storage in the computer 
than conventional elimination procedures. Nevertheless, a standard elimination sub- 
program was used. 
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Because of matrix-size limitations in the elimination subprogram used to solve 
equation (59), the technique of varying grid size (ref. 19) was used. This technique is 
illustrated in figure 8, where the grid spacing Ay is decreased toward zero. The at- 
tenuation at the optimum impedance is found by extrapolating the grid size to zero 
(dashed portion of the curve). As shown in figure 8, the extrapolated numerical values 
are in agreement with the analytical values calculated from the theory presented by Rice 
in reference 4. 

Based on the known computational times for the solutions presented in this report, 
we can now estimate the computer time involved in solving some realistic liner prob- 
lems. We will assume that an elimination subprogram will be available which is spe- 
cially tailored to the particular acoustic matrix; that is, no zero need be stored in the 
subprogram. Furthermore, we will assume that 1000 grid points will adequately de- 
scribe the liner geometry to be analyzed. For the two-dimensional example problem 
considered in this report, the solution would require approximately 0.4 minute of com- 
puter time. In the future, when nonuniform shear flow effects are incorporated into the 
governing flow equations, the computer time will increase by a factor of 3, to 1. 2 min- 
utes. Both time estimates are quite acceptable for liner design calculations. 

The numerical procedure is quite flexible and can be used to obtain details of the 
acoustic field structure, such as pressure, velocity, and intensity, at any position in the 
duct. Complications such as variable impedance or wall geometry are easily pro- 
grammed. 


CONCLUDING REMARKS 

Finite difference solutions are presented for sound propagation in a one -dimensional 
hard-wall duct and a two-dimensional soft -wall duct for zero Mach number. The results 
show the numerical procedure to be in agreement with the corresponding exact analytical 
results. Because the solution matrix for the acoustic flow field is not positive definite, 
conventional iteration techniques cannot be used to solve the difference equations. Before 
the finite difference formulation can be applied to more complicated problems which re- 
quire a large number of grid points, such as the sonic inlet, iteration techniques will 
have to be developed to overcome present grid-size limitations. 

The finite difference formulation is flexible and should be a powerful tool in the so- 
lution of more realistic studies of inlet and exhaust ducts of turbofan engines. The 
present formulation allows complete freedom in choosing the inlet pressure profile and 
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the complex impedance along both boundaries. The extension of the present formula- 
tion to both uniform flow and shear flow is straightforward . 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 10, 1973, 

501-04. 
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APPENDIX A 


SYMBOLS 


A 



a l’ a 2’ a 3 ,a 4 

C 

c 


c l> c 2’ c 3 

AdB 

E 

F 

f 

£ ( R j ) 

f(y) 

G(y) 

H 

I 


I’ 


coefficient matrix, eq. (58) 

submatrix, eq. (59) 

elements of matrix 

submatrix, eq. (59) 

speed of sound 

elements of matrix C 

decrease in decibels, eq. (44) 

acoustic power, eq. (43) 

matrix, eq. (58) 

frequency 

function of tL 

function of y 

arbitrary function of y 

channel height, fig. 1 

dimensionless acoustic intensity, eq. (41) 

dimensional acoustic intensity, eq. (39) 


v^r 


k 

L 

m 

n 

O 

P 

P' 

P 


complex propagation constant 
wave number, oo/c 
dimensional length of duct 

total number of grid rows (points in y -direction) 
total number of grid columns (points in x -direction) 
order of ‘ 

pressure column vector, eq. (58) 
dimensional pressure, P'(x’,y’,t) 
dimensionless pressure fluctuation, eq. (10) 



Pa 

p’ 

t 

u 

u’ 

X 

Ax 

x’ 

y 

Ay 

y' 

z 

K 

n 

K 

X 

p 

a 

9 

X 

u> 

Operators: 

v 2 

V 

Subscripts: 

e 

i J 
i 

t 


amplitude of pressure fluctuation at duct entrance 

dimensional pressure fluctuation 

time 

dimensionless acoustic particle velocity, eq. (22) 
dimensional acoustic particle velocity 
dimensionless axial coordinate, eq. (9) 
axial grid spacing 

dimensional axial coordinate, fig. 1 
dimensionless transverse coordinate, eq. (8) 
transverse grid spacing 
dimensional transverse coordinate, fig. 1 
acoustic impedance 

dimensionless specific acoustic impedance 

dimensionless frequency, eq. (12) 

acoustic conductance ratio 

wavelength 

density 

acoustic susceptance ratio 
dimensionless specific acoustic resistance 
dimensionless specific acoustic reactance 
circular frequency 

Laplacian 

gradient 

exit condition 

axial and transverse indexes, respectively, fig. 2 
mode index 

transverse y -direction 
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X 

axial position 

y 

transverse position 

Superscripts: 


! 

dimensional quantity 

* 

complex conjugate 

r ) 

vector quantity 

(i) 

real part 

(2) 

imaginary part 



APPENDIX B 


ONE-DIMENSIONAL SOUND PROPAGATION IN AN INFINITE DUCT 

Pressure Field 

For one -dimensional sound propagation in the x-direction, the dimensionless wave 
equation, equation (11) , reduces to 


£ 

dx 


(2777?) 2 P = 0 

2 


(Bl) 


In this case, the waves are uniform with no gradients in the y-direction. The general 
solution of equation (Bl) is 


p = Ae“ i27n7X + Be +i27n?x 


(B2) 


The boundary conditions for equation (B2) are as follows 


B = 0 (no reflecting waves) 


p = 1 at x = 0 


(B3) 

(B4) 


Applying conditions (B3) and (B4) to equation (B2) gives as a solution 


p - e 


■i2nr]x 


(B5) 


or 


p = cos 27T77X - i sin 2i jtjx 


(B6) 


Thus, in terms of the real and imaginary pressures in the numerical solutions 


p^ = cos 2 t7?7x 
( 2 ) 

p v ' = -sin 2nrjx 


(Bl) 

(B8) 
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Terminal Impedance 

The specific acoustic impedance £ at any point is defined by equation (25) and is 
given for the x -direction as 

r = -i27T77 JL (B9) 

3p\ 

0X 4 

Substituting equation (B5) into equation (B9) gives the impedance at any position x as 

c x = 1 (BIO) 

or the impedance at any position from equation (24) as 

Z = pc (Bll) 

In terms of the acoustic conductance and susceptance ratios, defined by equation (27), 

k = 1 (B12) 

a = 0 (B13) 

These results can be used in the numerical analysis to allow a finite length of duct 
to represent an infinite duct by using k = 1 and cr = 0. 

c c 


Intensity 

The dimensionless axial acoustic intensity is given in the body of this report as 


I = _L p (2) 9P^ . p (l) 9P (2) 
2717] 3x 3x 


(42) 


Substituting equations (B7) and (B8) into equation (42) and collecting terms gives 

I = 1 (B14) 

Thus, the axial acoustical intensity is invariant along the length of the duct. 
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APPENDIX C 


FINITE DIFFERENCE EQUATIONS 


The basic difference equations are determined by integrating the appropriate differ- 
ential equation over a unit cell, as shown in figure 9. The cell is enclosed by the dashed 
lines which are spaced midway between the grid lines (not shown). The grid lines would 
go directly through the grid points. In this particular problem, the governing equation 
is the wave equation as given by equation (14) or (15) in the body of the report. Thus, the 
finite difference integrals become 


f f V 2 p^ + (2 tt?7) V 1 ) 

J cell 


dA = 0 


vV 2 ) + mV 2 ) 


dA = 0 


cell 


By applying Green’s theorem for the plane region, these equations reduce to 



(Cl) 


(C2) 


(C3) 


(C4) 


where dA is a differential area element of the unit cell, ds is the differential length on 
the boundary of the unit cell, and n is the normal to the cell boundary. 

As usual, the difference formulation approximates the first derivatives as follows: 


OP 

0X 


p j+i.r p i.] 


i+ j 
2 


Ax 


1 + O Ax 2 
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0p 

3y 



iJt l 




Ay 



Wave Equation 

The wave equation in finite difference form, equation (47) in the body of the report, 
is developed by applying equation (C3) to the cell marked 1 in figure 9 . In evaluating the 
surface integrals, the gradients are assumed to be constant on each side of the cell and 
the interior pressure is assumed to be constant in the interior of the cell. Thus, equa- 
tion (C3) becomes 


pPI 


-p 


( 1 ) 


( 1 ) 


hi 


Ax 


Ay 


P- • i 

+ -±drl 


- P U) - P^ ■ 
h i Ax + 1+ ld 


Ay 


i “ p i i 
J Lu. Ay 


Ax 


pU) . P (1) 

+ .iliti hi Ax + (27I77) 2 P^ 1 ^. Ax Ay = 

Ay 1>] 


(C5) 


Multiplying all terms by -Ax/Ay and collecting terms gives 


-p! 1 ) . 


P i+1 , j 


Ax 

Ay, 


p^ 1} . t + p< 1} . t 

^1, 3-1 


2 + 2 (— f - (2 ttt 7 Ax) 2 
Ay/ 


- 0 


(C6) 


which is equation (45) in the body of the report. Using the same procedure, the differ - 

( 2 ) 

ence equation for the p v ' component can also be found. 

Finally, let us consider the difference equation which applies in cell 2, which is 
adjacent to the upper boundary in figure 9. For this unit cell, equation (C3) can be ex- 
pressed as 
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_(l> - D ( " P (1> , - P U) P ( H - P^' a 

p i-l,m p i.m Ay + p i,m-l p i,m ^ t p i+l,m i,m Ay 

Ax 2 


Ax 


Ay 


dp 


( 1 ) 


3n 


wall 

y=i 


Ax + (2 ttt?) 2 p^ Ax = 0 
1 > m 2 


(C7) 


The normal derivative at the wall, 3p/3y is given by the known impedance boundary con- 
dition given by equation (29) in the body of the report, which is repeated here for con- 
venience 


= -2 Jp<\ - p®« t ‘ 
3y L t tj 


(29) 


In difference form, 




( 1 ) 


3y 


= -27777 


y=i 




m 


(C8) 


Substituting equation (C8) into equation (C7) and collecting terms yields 


P 


( 1 ) 


i, m 


Ay 

Ax 


1 + (~| + 2nrj AylcrJ 


(2777 ) Ay Y 


l, m 


-p(l> i/aAm 

i,m-l 2 \Ax/ 1-1 ’ 


m 


- ; i ( ?i , m - 2 «> Ay(«t) p® = 0 <C9> 

1 \Ax/ ’ ' m ’ 


In a similar manner, equations (C3) and (C4) are applied to cells 2, 3, and 4 to 
yield equations (50) to (56) in the body of this report. 
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APPENDIX D 


APPLICATION OF DIFFERENCE EQUATIONS TO SOFT-WALL DUCT 

In this appendix, we develop some of the matrix coefficients which make up the 
matrix given by equation (60) in the body of this report. This matrix results from the 
difference formulation for sound propagation in a straight soft-wall duct for zero Mach 
number . 


General Difference Equations at Lower Wall 


Before looking at the specific form of equation (60) , we will write the general form 
of the difference equations at the lower wall, since the coefficients in the first two rows 
of equation (60) result from applying the difference equations to points l and 2 along the 
lower wall shown in figure 3. At the lower wall, equations (48) and (49) can be written 
as 


p) 1 ! 

I 1 


1 + 


M 2 * 2»„ AyU - . p<« . I (*lf p<» 

Ax/ ' V Vi, 1 2 1)Z 2 \A x/ 11)1 


j(^) 2 pw,i - 2 ”" = ° < di > 


p 


( 2 ) 

i, 1 


1 + 


(— ^ + 27777 Ay (a,'] 
\Ax/ V Vi I 


(2nij Ay)‘ 
2 


(2) 1 / AyV (2) 




I /Ay' 

2 \Ax 


p (2) 

p i+l , 1 


2nr] Ay 



= 0 


(D2) 


In this example, Ay = Ax; therefore, the ratio Ax/ Ay = 1; thus, these equations reduce 
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p 


( 1 ) 

i. 1 


2 + 27777 Ay 


(°‘), r 


(2777? Ay)‘ 


0 (« .1 

1 - 2 2 


D (l) + B (D 
Pi-1,1 + p i+l,l 


- 2777] Ay^c^ ^p| 2 ^ = 0 (D3) 


3 (2) 

i.l 


2 + 2777] Ay 


R,r 


(2777? Ay)‘ 
2 


p(2) - I 
i,2 2 




+ 27777 Ay 


R/m ' 0 


(D4) 


For simplicity, we now define 


aj = 2 + 2777] Ay^ 


(2777] Ay)' 
2 


(D5) 


and 


c x = 2777] Ay bcj 




(D6) 


Substituting equations (D5) and (D6) into equations (D3) and (D4) yields 


a x> (l) - o (1) 1 

a l P i,l P i , 2 - - 


} (D + p (i) 

5 i-l, 1 + p i+l, 1 


- c p^ = 0 
c l p i, 1 u 


(D7) 


a n (2) - d (2) 1 

a l p i,l p i,2 " 


p (2) + p (2) 

p i-l , 1 + p i+l , 1 


+ c p(D 

1 i 1 = 0 


(D8) 


In a similar manner, the difference equations at the lower corner take on the form 


a p(D _ D (l) _ p(D _ c d ( 2) _ o 

a 2 p n, 1 p n-l , 1 p n, 2 c 2 p n,l~° 


a p {2) - p( 2 ) - p (2) + c n (1) - 0 

a 2 P n, 1 p n-l , 1 P n,2 C 2 p n,l _U 


(D9) 

(DIO) 
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where 


a 2 = 2 + 27TT7 Ay(a t ^ ^ - 


c 2 = 2rq Ay 


1 + K. 


(2irr] Ay Y 
2 


n, 1 


(DU) 

(D12) 


Application to Equation ( 60 ) 


We shall now apply these equations and the other equations derived elsewhere in this 
report to the derivation of equation (60) . 

Point 1 . - The double-subscript notation will be transformed into the single- 
subscript notation used in figure 3 for this example. Applying equation (D7) to point 1 
yields 




P (1) 

p 0 




(D13) 


where a^ and Cj are defined by equations (D5) and (D6) . However, p Q represents an 
initial -condition grid point, which in this example has a value of 1. Thus, equation 
(D13) becomes 


aiP 


( 1 ) 


1*1 




vP-5 


(D14) 


or in matrix form 
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(D15) 


which represents the first row in equation (60). 

The imaginary pressure equation at point 1 results from applying equation (D8) 

to point 1 . This yields 


In this case, however, 
comes 


aiPf - Pi 2 ’ - \ 


pS 2) - 4 2) 


+ CjP^ 1 ) = 0 


(D16) 


Pq 2 ^ represents a zero initial condition; thus, equation (D16) be- 


,(D 


iPi + a 1?i 


(2) 


M 2) 

2 ^ 


p^ 2) = 0 


(D17) 


which represents the 10th row of equation (60). The matrix form of equation (D17) is 
similar to equation (D15). 
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Point 2. 

- At point 2, equations (D7) and (D8) become 



-ip) 1 ' t ajp) 1 ' - ip) 1 ' - p) 1 ' - c lP ® = 0 

(D18) 

and 

-ipf'.a^'-ipf-pf.c^'^O 

(D19) 

which represents the second and 11th rows of equation (60). 
Point 3. - At point 3, equations (D9) and (DIO) become 



-P^ + a 2P3 J) - p 6 } - c 2 p 3 2) = 0 

(D20) 

and 

" p 2 2) + a 2 p 3 2) - p 6 2> + c 2 p 3 1} = 0 

(D21) 

which represent the third and 12th rows of equation (60). 

Point 4. - At point 4, equations (45) and (46) in the body of this report become 



-P^ - P ( 5 X) - Pj 1 ^ - P ^ + 4 - (27rq Ax) 2 p^ - 0 

(D22) 

and 




~ p 0 2) - p 5 2) - p l 2) - p 7 2) + 4 " ^ p 4 2) = 0 

(D23) 

But from the initial condition, = 1 while p^ = 0; thus equations (D22) and (D23) 

become 


-p<" + a 3 p< 1 ' - p) 1 ' - p' 1 ' . 1 

(D24) 


-pf.a 3 pf -P< 2 '-P® = 0 

(D25) 
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where 


a^ = 4 - (27r?7 Ax)^ 


(D 26 ) 


Equations (D24) and (D25) represent the fourth and 13th lines in equation (60). 

Point 5 . - The difference equations that apply at point 5 come directly from equa- 
tions (45) and (46) as 


and 






(D27) 




+ a 3>4 2) 




(D28) 


These equations represent the fifth and 14th lines in equation (60). In an actual acoustic 
problem, the bulk of the difference equations in the central flow field will be of this form. 
Point 6. - At point 6, equations (52) and (53) in the body of this report become 


and 


-jpSs 1 ’ - 1>5 ’ + Ve’ - 1 Pg 1 ’ - 


c 3 p® = 0 


(D29) 


1 

2 




+ a 4 p 6 2 ' 


ip® 

2 y 


SD- 

'3 P 6 


+ c 0 Pp - ' = 0 


(D30) 


where 


a _ 2 _ (27777 Ax) 2 


(D31) 


and 


Cg = 27777 Ax (D32) 

Equations (D29) and (D30) represent the sixth and 15th lines in equation (60). 

The equations for points 7, 8, and 9 are similar to the equations at points 1, 2, 
and 3. These equations will not be presented here. 
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APPENDIX E 


CLOSURE PROBLEM 


The close agreement between the numerical theory and the corresponding exact 
analytical results substantiates the utility of the pc exit impedance assumption made 
earlier in this report. However, we shall now establish criteria for which the numerical 
and analytical results will be in exact agreement. To accomplish this task, we must 
first look at the analytical theory. 

The analytical solution (ref. 7, e.g.) assumes that the solutions to the wave equa- 
tion are separable and can be expressed in terms of an infinite number of modes. For 
only forward-going waves (no reflections, i.e. , an infinite duct) the expressions for 
pressure and velocity are of the form 


N 

Z -ikKkx 
G.(y)e J 

j=l 


(El) 


where K. is a complex propagation coefficient for each mode. 

The propagation constants K. are functions of an eigenvalue p... The eigenvalues 
\ * J ] 

p.j are determined by the wall impedance condition 


£ t = 27777 




M 


w 


(E2) 


A solution to equation (E2) yields an infinite set of eigenvalues ju^, each identifiable with 
a mode. Also, as a direct consequence of equation (E2), it follows that for each mode 



Conversely, it follows that at some finite distance L down the infinite duct, which 
we label e for exit plane, the ratio of pressure to axial velocity is different for each 
mode : 
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(E4) 



Equation (E4) follows as a consequence of the fact that this ratio depends on the eigen- 
values (through K.) which were previously determined by equation (E2): 



(E5) 


Thus, each ratio of pressure to axial velocity at the exit position takes on its own unique 
value depending on the values of the propagation constant K. . 

If we now define the axial impedance in the usual manner, the impedance condition 
for no reflection at x = L becomes 


P 1 p 2 Pi . . 

£ 1 = 2W , f „ = 2tt77 , . . . , £ • = 2tt7? —L, . . . (E6) 

M ”x2 u xj 

That is,' the infinite duct with no reflections could now be replaced by a finite duct of 
length L with each mode having its own respective exit impedance as given in equa- 
tion (E6). Clearly,' as a consequence of equations (E5) and (E6), the exit impedance con- 
ditions for no reflection are such that 

S el ^ S e2 ^ S e3 ^ ^ £ ej ^ ^ 

with each mode having its own unique exit impedance for no reflection. 

In the numerical treatment, the boundary condition used at the duct exit x = L was 

£ = 27777 -£•= 1 (E8) 

" ”x 

where the numerical value of 1 indicates the condition for the case in which a plane one- 
dimensional pressure wave would not be reflected at the tube exit. 

In comparing the exact analytical exit impedance conditions required for no reflec- 
tion, equation (E6), to the numerical exit impedance condition, equation (E8), we are led 
to the conclusion that reflections are occurring at the duct exit in the numerical calcula- 
tions. Choosing equal to 1 in the numerical calculation cannot satisfy all the distinct 
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impedance values for each mode as given in equation (E6). Clearly, the analytical and 
numerical solutions for the attenuation in a duct of a given length will not be identical. 

If we overlook roundoff and truncation errors in the numerical calculations, the at- 
tenuation calculated by the numerical technique will be greater than that calculated by 
the exact analytical technique because of reflections which occur at the exit; that is, 

numerical ^ analytical 

0<x<L 0<x<L 


The reflected energy appears to have been absorbed by the soft wall in the numerical 
calculation. 

Fortunately, in many practical problems, most of the higher order modes decay out 
and only the lowest order mode appears at the exit. A comparison to the analytical re- 
sults indicates that in general the exit impedance associated with the lowest order mode 
is close to pc. Nevertheless, it is important to establish a general numerical proce- 
dure which will converge to the correct answer. 

Consider the case where an additional length of absorbing liner AL is added to the 
original liner, so that the liner now has a length L + AL. If AL is chosen sufficiently 
long, the reflected energy at the duct exit, x = L + AL, will be absorbed before this 
energy can reach the previous exit position at x = L. Therefore, for soft -wall ducts 


limit AdB 
AL — 00 


numerical 

0^x<L 


) analytical 
0<x<L 


(E10) 


Practically, the liner length need only be increased by some minimum length AL until 
the attenuation for the liner length between x = 0 and x = L remains constant to a 
given percentage. 


39 



REFERENCES 


1. Morse, Philip M, : The Transmission of Sound Inside Pipes. J. Acoust. Soc. Am., 

vol. 11, no. 2, Oct. 1939, pp. 205-210. 

2. Cremer, Lothar: Theorie der Luftschall-Dampfung in Rechteekkanal mit 

Schluckender Wand und das sich dabei Ergebende hOchste D£mpfungsmass. 
Acustica, vol. 3, no. 2, 1953, pp. 249-263. 

3. Rice, Edward J. : Attenuation of Sound in Softwalled Circular Ducts. Presented at 

the AFOSR-UTIAS Symposium on Aerodynamic Noise, Toronto, May 20-21, 1968. 

4. Rice, Edward J. : Propagation of Waves in an Acoustically Lined Duct with a Mean 

Flow. Basic Aerodynamic Noise Research. NASA SP-207, 1969, pp. 345-355. 

5. Eversman, Walter: The Effect of Mach Number on the Tuning of an Acoustic Lining 

in a Flow Duct. J. Acoust. Soc. Am., vol. 48, no. 2, pt. 1, 1970, pp. 425-428. 

6. Tack, D. H. ; and Lambert, R. F. : Influence of Shear Flow on Sound Attenuation 

in a Lined Duct. J. Acoust. Soc. Am. , vol. 38, no. 4, Oct. 1965, pp. 655-666. 

7. Pridmore- Brown, D. C. : Sound Propagation in a Fluid Flowing Through an Atten- 

uation Duct. J. Fluid Mech. , vol. 4, pt. 4, Aug. 1958, pp. 393-406. 

8. Mungur, P. ; and Gladwell, G. M. L. : Acoustic Wave Propagation in a Sheared 

Fluid Contained in a Duct. J. Sound Vib. , vol. 9, no. 1, Jan. 1969, pp. 28-48. 

9. Mungur, P. ; and Plumblee, H. E. : Propagation and Attenuation of Sound in a 

Soft- Walled Annular Duct Containing a Shear Flow. Basic Aerodynamic Noise 
Research. NASA SP-207, 1969, pp. 305-327. 

10. Hersh, A. S. ; and Catton, I.: Effect of Shear Flow on Sound Propagation in 

Rectangular Ducts. J. Acoust. Soc. Am., vol. 50, no. 3, pt. 2, 1971, pp. 992- 
1003. 

11. Eversman, Walter: Effect of Boundary Layer on the Transmission and Attenuation 

of Sound in an Acoustically Treated Circular Duct. J. Acoust. Soc. Am. , vol. 49, 
no. 5, pt. 1, 1971, pp. 1372-1380. 

12. Mariano, S. : Effect of Wall Shear Layers on the Sound Attenuation in Acoustically 

Lined Rectangular Ducts. J. Sound Vib. , vol. 19, no. 3, Dec. 8, 1971, pp. 261- 
275. 

13. Savkar, S. D. : Propagation of Sound in Ducts with Shear Flow. J. Sound Vib. , 

vol. 19, no. 3, Dec. 8, 1971, pp. 355-372. 

14. Shankar, P. M. : On Acoustic Refraction by Duct Shear Layers. J. Fluid Mech. , 

vol. 47, pt. 1, May 14, 1971, pp. 81-91. 


40 



15. Ko, S. -H. : Sound Attenuation in Acoustically Lined Circular Ducts in the Presence 

of Uniform Flow and Shear Flow. J. Sound Vib. , vol. 22, no. 2, May 22, 1972, 
pp. 193-210. 

16. Shankar, P. N. : Sound Propagation in Duct Shear Layers. J. Sound Vib. , vol. 22, 

no. 2, May 22, 1972, pp. 221-232. 

17. Shankar, P. N. : Acoustic Refraction and Attenuation in Cylindrical and Annular 

Ducts. J. Sound Vib. , vol. 22, no. 2, May 22, 1972, pp. 233-246. 

18. Varga, RichardS.: Matrix Iterative Analysis. Prentice- Hall, Inc., 1962. 

19. Moon, Parry; and Spencer, Domina E. : Field Theory for Engineers. D. Van 

Nostrand Co. , Inc., 1961. 


41 




x' ■ L 


Figure 1. - Coordinate system for two-dimensional duct. 




grid points 


Figure 2. - Grid-point representation of two-dimensional flow duct. 
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Figure 3. - Illustrative example of soft -wall duct with total number of grid rows m of 3 and 
total number of grid columns n of 3. 



Figure 4. - One-dimensional sound propagation with no reflections. 
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Dimensionless axial intensity at duct exit, 


Theory (eq. (B14)> 







Transverse grid spacing, Ay = l/(m - 1) 

Figure 8. - Effect of spacing on attenuation at opti- 
mum impedance. 
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Figure 9. - Integration cells for establishing difference equations. 
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